Spin flop in one-dimensional quantum antiferromagnets 
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The field-induced transition corresponding to the spin flop in the quantum antiferromagnetic 
chains in the presence of the Ising-like single-ion anisotropy D (< 0) is studied by the finite-cluster 
analysis. It is found that for S > § a first-order metamagnetic transition occurs even in one 
dimension for small negative D except for the Haldane phase, while there exist two second-order 
transitions and an incommensurate phase appears between them for intermediate D , as already 
found for S = 1. We also discuss on some experiments of CsNiCl3 related with the present work. 

PACS numbers: 75.10.Jm, 75.40Mg, 75.50.Gg, 75.40.Cx 



Magnetization measurements of antiferromagnets de- 
tect various macroscopic quantum effects. One of inter- 
esting examples appear in the spin flopping (SF) process 
of low-dimensional materials. SF is one of field- induced 
metamagnetic phase transitions of anisotropic Heisen- 
berg antiferromagnets. [jj] In terms of the classical spin 
systems, the transition brings about an abrupt change in 
the direction of the Neel order from parallel to perpen- 
dicular with respect to the easy axis under the applied 
field H along the axis and it is a first-order transition 
with a jump in the magnetization to as a function of H. 
In quantum systems the jump generally shrinks due to 
quantum fluctuation and such an effect becomes larger in 
lower dimensions and for smaller S. For S = ^ antifer- 
romagnets with the Ising-like anisotropic exchange cou- 
pling, the magnetization jump was revealed to survive 
down to two dimension by the recent numerical analy- 
sis |g], while the exact solution || by the Bethe ansatz 
suggested that it changes into a second-order transition 
descibed by the asymptotic behavior 
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where 8 = 2 in one dimension (ID). Thus the large quan- 
tum fluctuation in ID seems to change the order of the 
transition in the system. 

In view of experimental realization, one of the most 
important anisotropy is the single-ion anisotropy (SIA) 
described by D J^ji^j) 2 which is relevant only for S > 1. 
In the classical limit in any dimensions the negative D 
makes the Neel order oriented along z-axis and the first- 
order SF transition is induced by the applied field parallel 
to the axis. The recent numerical analysis B suggested 
that instead of SF two successive second-order transitions 
occur in ID S = 1 antiferromagnet for small negative D 
except for the Haldane phase. Q This is because the 
magnetic phase of the system consists of the two gap- 
less Luttinger liquid phases with different elementary ex- 
citations created by (Sj~) 2 and Sj~, respectively. The 
first transition, described by the form (||) with 6 = 2, 
occurs from the Neel state to the first massless phase 
characterized by the dominant spin correlation function 
{SqS*) ~ cos{2k F r)r~ r '% where 2k F = 2wm, and the 
next transition with (5 = 1 leads to the second massless 



phase where (Sq S~) ~ (— l) r r~ n . The behavior of the 
spin correlation function suggests that the first massless 
phase, which lies between the two transitions, has no cor- 
respondence in the classical systems. It implies that the 
appearance of the phase is an intrinsic quantum effect. 
Is such a strongly quantized behavior general in ID, or 
peculiar to S — 1? The question encourages us to investi- 
gate ID larger-,!? systems. In this paper, using the finite- 
cluster analysis, we study the magnetization process of 
the ID S — | and 2 antiferromagnets in the presence of 
the Ising-like SIA, in comparison with the case of S = 1. 
The ID Heisenberg antiferromagnet with SIA under 
magnetic field parallel to the easy axis of the Neel order 
is described by the Hamiltonian 

H = Ho + Hz , 

H =J2 S rS 1+ i + Dj2( S j) 2 i (2) 

i i 

Hz= — H 2_^ Sjt 



where S^ = S(S+1) and the periodic boundary condition 
are applied. We restrict us on the Ising-like anisotropy 
D < to investigate the transition corresponding to SF. 
For L-site systems, the lowest energy of Ho in the sub- 
space where the eigenvalue of V S? is M (the bulk mag- 
netization is to = M/L) and the momentum is k, is de- 
noted as Ek(L,M). In addition we define E(L,M) as 
the lowest one among Ek(L,M)'s. Using Lanczos algo- 
rithm, we calculated Ek(L,M) for even-site systems up 
to L = 20, 14 and 12 for S = 1, § and 2, respectively. 

The magnetic ground state (GS) for S > 1 is divided 
into two massless phases B in the phase diagram on the 
rn/S-D (D < 0) plane as well as S — 1. The argument 
on the two phases for S = 1 developed in the Ref.[4] is 
easily generalized for arbitrary S as follows: (i)Large neg- 
ative D phase; the gapless excitation is created by (S^~) 2 
and the soft mode has the momentum 2kp = (1 — ^§)tt. 
(ii)Small negative D phase; the S^ excitation is gapless 
and the soft mode has 2k' F = 2mSir. The dominant spin 



correlation functions are (i)(SqS^) ~ cos(2fci?r)r Tlz and 
(ii)(SqS~) ~ (— 1)'V~'', respectively. Thus in quasi-lD 



the coherent interchain interaction is expected to result 
in the long-range (i)incommensurate SDW along the z- 
axis and (ii)Neel order in the xy-plane (called the canted 
Neel order in some texts), respectively. In this paper we 
call the two massless phases (i)'z-SDW' and (ii)'a;?/-Neel' 
phases to emphasize the long-range order expected to be 
observed in real quasi-lD systems, whereas pure ID sys- 
tems have some power-low decays in the associated spin 
correlation functions. 

Consider the two energy gaps defined as 



Ai = E(L, M + 1) + E(L, M • 
A 2kF =E 2kF (L,M)-E(L,M). 



1)-2E(L,M), 



(3) 



Ai is the sum of the gaps corresponding to SM = 1 and 
— 1 for the total Hamiltonian 7i. Since Ai (A2k F ) is open 
(gapless) in the z-SDW phase, while gapless (open) in the 
xy-Neel phase, the crossing point of the two gaps for fixed 
to is a good estimation of the phase boundary. Since only 
states with M — 2Sn (n = 0, 1, 2, • • •) contributes to the 
magnetization process in the z-SDW phase, our analysis 
is focused on such states. The crossing points of the fi- 
nite clusters with various magnetizations are plotted on 
the m/S-D plane as open symbols in Figs. 1 (a)5 = | 
and (b)5 = 2. Little size dependence of the curves justi- 
fies that the fitted solid lines are the phase boundary of 
the bulk systems. D c q denotes the boundary between the 
Neel ordered and disordered phases in the nonmagnetic 
GS. It corresponds to the boundary of the Haldane phase 
for integer 5, while D c q — for half-odd integer 5. Us- 
ing the phenomenological renormalization pj , D c q of the 
5 = 2 system is estimated as D c0 = -0.001 ±0.001. The 
phase boundary at to/5 = 1 denotes D C 2- It was deter- 
mined as the crossing point of E(L, SL) — E(L, SL — 1) 
and E(L, SL) — E(L, SL — 25) which is almost indepen- 
dent of L. The results are D C 2 = —1.31 and —1.25 
for 5 = | and 2, respectively. The m/S-D phase di- 
agrams suggest that besides the two ordinary critical 
fields H c i and H C 2, which are the starting and satura- 
tion points of the magnetization, respectively, the mag- 
netization process has an intermediate critical field H c ^ 
(H c i < H c3 < H C 2) corresponding to the boundary of 
the two massless phases for D c0 > D > D c2 . Accord- 
ing to the argument for 5 = 1, the phase transitions at 
the three critical points are all second-order and have the 
asymptotic forms 



(H - H cl 



.1/2 



m - m c3 ~ \H - H, 



<-;', 



(H C 2 Ufl\ 



(4) 



(5) 



(6) 



respectively, where to C 3 is the magnetization correspond- 
ing to the boundary between the two massless phases. In 
the classical limit the first-order SF transition occurs and 
to jumps from to [— D/(D + 2)] 1 / 2 5 (long-dashed lines 



in Figs. 1 (a) and (b)) for > D > -1, while to 5 
for D < — 1. The phase diagrams including only the two 
massless phases might lead to the conclusion that even 
for larger 5, as well as 5 — 1, two second-order transi- 
tions occur instead of SF for D c q > D > D c2 and the 
phase boundary between the two massless phases will go 
to the classical line as 5 increases towards oo. 



D,,=-1.31 



1.0 



0.6 



0.4 



0.2 



0.0 



(a)5=3/2 I 

o L=14 4 
■ □ L=12 % 


\ 


a L=10 \ 


\ 


• £"(m)=0 % 


\ 




xy-Nete| 


z-SDW 






\ e"(m)<T\^V 



-2.0 



1.5 


-1.0 


-0.5 


0.0 


D 


t 




T 




D„=-0.9 




Z>,„=0 



D„=-1.25 




1.5 


-1. 


D 


T 




D =-1.1 



0.0 

t 

D co =-0.00l 

FIG. 1. m/S-D phase diagrams for (a)5 = § and 

(h)S = 2. Open symbols are the crossing point of the two 
gaps Ai and A2fc F for finite systems. The fitted solid curves 
are the boundaries between the z-SDW and rry-Neel phases. 
Solid circles are the points with e"(m) = and the curvature 
e"(m) is negative in the regions surrounded them. (Dashed 
curves are guides for the eyes.) 



In the following analysis, however, it will be found that 
the first-order SF transition can occur for 5 > |. The 
lowest eigenvalue of the Hamiltonian Ho for < to < 5 
per site in the thermodynamic limit, denoted as e(m), 



can be estimated using the form 



jE(L,M) 



e(m) 



Cim)^, 



(J) 



where the second term is the size correction. As the GS 
magnetization curve is derived from H = e'(rn), a nec- 
essary condition of the continuity of the curve at to is 
e"(m) > 0. Since only M = 2Sn is available in the z- 
SDW phase, we investigate the condition using the form 



R (L, M) = 
L [E(L, M 

~ e (to) 



25) + E{L, M - 25) 



2E(L,M)]/(2S) 2 
(8) 



derived from (M) 

3 

2 



#(£,25) is plotted vs D for 5 = 1, 
and 2 in Figs. 2 (a), (b) and (c), respectively. They 
indicate the existence of the region with R(L, 25) < 
for 5 = | and 2, and it tends to extend towards D = 
with increasing L, in contrast with the case of 5 = 1 
where R(L, 25) is always positive. Neglecting the size 
correction in (||) and taking R(L, M) for the largest L 
available for each to as e"(m), the points with e"(m) = 
are plotted as solid circles in Figs. 1 (a) and (b). In 
the area surrounded by them e"(m) is negative. The fit- 
ted dashed lines are guides for the eyes drawn assuming 
that the negative e"(m) region continues to D c q. Only 
for 5=1, e"(m) > is satisfied everywhere on the m-D 
plane within the same analysis up to L — 20. For 5 = § 
and 2 the existence of the area with e"(m) < implies 
that the first-order SF transition occurs and to jumps 
over the region. Figs. 1 (a) and (b) suggest that SF can 
appear only in the small negative D region. This is be- 
cause the large negative D strongly quantizes the system 
and makes it equivalent to the 5 = \ XX Z model which 
is proved to have no first-order transition. Obviously the 
negative e"(m) region is larger for 5 = | than 5 = 2, 
which is consistent with the assumption that smaller- 5 
systems have larger quantum fluctuation suppressing the 
magnetization jump. The large negative D border of the 
SF region, denoted as D cl , is roughly estimated from the 
boundary of the negative e"(m) region for the finite sys- 
tems as shown in Figs. 1 (a) and (b). The results are 
D c \ ~ —0.9 and —1.1 for 5 = | and 2, respectively. The 
size correction neglected on the estimation of D c \ is not 
supposed to be so large as to change the following qual- 
itative argument, because Figs. 2 (b) and (c) suggest 
that the left border of the region with R(L, M) < has 
no significant L dependence in comparison with the right 
border. For D c q > D > D c \ the magnetization process 
is expected to have a jump from to = to the xy-Neel 
phase. Figs. 1 (a) and (b) might suggest the possibility 
of the jump from nonzero to near Z? c i, but it is difficult 
to determine the point of the jump within the present 
small cluster analysis because the first order transition 
is ruled by the Maxwell construction, which cannot be 
applied until a more precise curve of e(m) is given. 
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FIG. 2. R{L, 25) vs D of the finite systems for (a)5 = 1, 
(b)5 = | and (c)5 = 2. R(L,2S) vs 1/L with fixed D; 
(d)5 = 1 D = -0.8, (e)5 = § D = -0.6 and (f)5 = 2 
D = —0.6. These values of D are around the minima (local 
minimum for S = 1) of the largest-i curves in (a), (b) and 
(c). Dashed lines are based on the least square fitting. 



To convince of the existence of SF for 5 = | and 2 
in the thermodynamic limit, we investigate the behavior 
of e{m) around the limit to —* 0+. Assuming that the 
size correction of e(m) is at most 0(1/L 2 ) even in the 
massive phase (to = 0), the size dependence of R(L,2S) 
is expanded around to = as 



R(L, 25) ~ e"(0) + (25)e"'(0)i + O(-l). (9) 



Around the minimum (local minimum for 5 = 1) of 
the R(L, 2S)-D curve in Figs. 2 (a-c) for the largest 
L (D = -0.8 for 5 = 1 and D = -0.6 for 5 = § and 
2), the size dependence of R(L, 25) is investigated by the 
plot of R(L,2S) vs 1/L in Figs. 2 (d)5 = 1, (e)5 = § 
and (f)5 = 2, where dashes lines are the results from the 
least square fitting. Considering the form (H), the plots 
for 5 = § and 2 both suggest e"(0) < and e"'(0) > 
which is a sufficient condition of the appearance of the 
magnetization jump from to = at some critical field in 
the thermodynamic limit. On the other hand, the plot for 
5=1 implies e"(0) = and e"'(0) > which is consis- 
tent with the second-order transition at H c \ described by 
the asymptotic form (Q). These results also support the 
existence of SF for 5 > | in the thermodynamic limit. 
Since the quantum fluctuation becomes smaller for larger 
5, the SF region should become larger with increasing 5 
and in the limit 5 — » oo the region would occupy all the 
region at the left hand side of the classical SF line in the 
m/S-D phase diagrams in Figs. 1 (a) and (b). 
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FIG. 3. Schematic magnetization curves of four possible 
cases. Solid and dashed curves signify the xy-Nee\ and z-SDW 
phases, respectively, except for massive regions. 



Based on the m/S-D phase diagrams in Figs. 1 (a) 
and (b), we summarize all possible cases of the magne- 
tization process of ID Heisenberg antiferromagnet with 
arbitrary S (> \) in the presence of the Ising-like SIA. 
In the following discussion the critical fields H c \ and H c2 
always stand for the starting and saturation points of 
the magnetization. The schematic magnetization curves 
are shown in Figs. 3 (a), (b), (c) and (d), and their 
features are listed as follows: (&)D > D c q (only for inte- 
ger S); the second-order transition with the asymptotic 
form (3) from the Haldane to xy-Nee\ phases occurs at 
H cl . (b)D cQ > D > D c i, the first-order SF transition oc- 
curs at H c i from the Neel to xy-Neel phase. The present 
analysis leads to the absence of the case only for S = 1. 
(c)D c i > D > D c2 ; the second-order transition from the 
Neel to z-SDW phases at H c \ described by the form (Q) 
is followed by another second-order transition to the xy- 
Neel phase at H c ^ with the form (H). (d)D < D c2 \ the 
second-order transition from the Neel to z-SDW phases 
and the latter phase continues until the saturation. The 
last case (d) in the large negative D region is equivalent 
to the Ising-like S — \ XXZ chain. In any case the 
transition at the saturation point H c2 is a second-order 
one described by the critical behavior (g). As S increases 
towards the classical limit, the boundaries D c \ and D c2 
will go to the same limit 1 and the case (c) will disappear 
in the limit. For D < D c2 H c \ and H c2 will become close 
to each other and coincide in the limit 5* — ► oo. 

The most interesting case is (c) in Fig. 3 where the 
two second-order transitions take place and the z-SDW 
phase appears between them. The most realistic system 
where such a magnetization process can be realized is 
the 5 = 1 chain. One of good candidates to observe the 
process is expected to be CsNiChj, a quasi- ID 5 = 1 



antiferromagnet with the Ising-like SIA, which has the 
Neel order at low temperatures. M The high-field mag- 
netization measurements |l(],[ll| on CSMCI3 indicated a 
SF-likc transition. However, the experimental magne- 
tization curve also looks like a second-order transition 
described by the critical behavior (Q), or two successive 
transitions accompanied by an intermediate phase. |12] 
In addition the NMR measurement |13| on CsNiCi3 indi- 
cated two spin structures oriented in different directions 
around the transition. It may imply that another tran- 
sition from the z-SDW to xy-Neel phases occurs soon 
after the first one. It would be difficult to detect the 
latter transition at H c ^ in the measured magnetization 
curve, because the transition brings about no significant 
anomalous behavior in the curve. || As well as NMR, the 
ESR or neutron scattering would be a better method to 
detect the change in the direction of the dominant spin 
correlation or the long-range order from parallel to per- 
pendicular with respect to H on the transition at H c3 . It 
would be interesting to observe the z-SDW phase which 
is an essential quantum effect. 

Since the quantum effect is generally larger for smaller 
systems, the first-order SF transition is more difficult to 
occur in finite clusters, treated in the present study, than 
in the infinite system. Therefore the conclusion that SF 
due to SIA can occur for S > | in ID is expected to be 
valid even in the thermodynamic limit. In contrast, the 
absence of SF for S = 1 should be checked by some other 
methods like the Monte Carlo or DMRG Q to treat 
larger chains. It would be also interesting to perform 
a magnetization measurement on some ID larger-S* ma- 
terials with negative D (\D\ < J), where the first-order 
transition is expected to occur, and to compare the curve 
with that of the S = 1 chain CsNiCl 3 . 

In summary the finite-cluster analysis on ID 5 = | 
and 2 antiferromagnets suggested that the first-order SF 
transition can occur for S > | in the presence of the 
Ising-like SIA. We also discussed on all possible magne- 
tization processes for D < including the case when two 
successive second-order transitions appear instead of SF. 
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